Quantification of group chasing and escaping process 
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We study a simple group chase and escape model by introducing new parameters with which 
configurations of chasing and escaping in groups are classified into three characteristic patterns. In 
particular, the parameters distinguish two essential configurations: a one- directional formation of 
chasers and escapees, and an escapee surrounded by chasers. In addition, pincer movements and 
aggregating processes of chasers and escapees are also quantified. Appearance of these configurations 
highlights efficiency of hunting during chasing and escaping. 

PACS numbers: 47.54.-r,05.65.+b,87.18.Gh,87.18.Hf 



Collective phenomena, such as school of fish and for- 
mation flight of birds, have been attracting interests for 
a long time. Recently, problems on various self-propelled 
particles have been extensively researched in the context 
of physics0,i[l3l. Previous works mainly focused on 
spontaneous group-formation in experiments and math- 
ematical models. Indeed, group motions often become 
crucial for survival in nature. For example, collective 
motions in hunting have advantages to chase a prey for 
a predator as well as to escape from a predator for a 
prey. Hunting is a widespread phenomenon observed 
from micro-organisms such as a neutrophil capturing 
viruses [l[ to macroscopic entities in animal societies 
and r oboticsd Q. Preys and predators adopt various 
patterns of collective motions as a result of competition 
for survival. A way to select an optimal strategy for each 
group is an interesting point. In a hunting motion, basic 
actions are chases and escapes. A problem of one-on-one 
chase-and-escape also has attracted mathematicians with 
its complex trajectories even in simple setups [3, 3- How- 
ever, group hunting could not be explained just by sum- 
ming up of individual performances of prey or predator 
and can be quite more complex. For example, coopera- 
tion within predator or prey group can give more diverse 
strategies for chasing or escaping. A pincer movement 
is one typical example. The emergence of collective mo- 
tions in group hunting has also been attracting interests, 
however quantitative study has not been done. 

One example of capturing and escaping patterns is 
shown in Fig. [T] which is actually a result of numerical 
simulation described later. Chasers surround aggregate 
of escapees to capture them. On the other hand, a mass 
of escapees gets away from chasers while they chase "de- 
coy" escapees. These patterns of group motions seem 
universal behaviors in hunting. If these collective mo- 
tions have an advantage rather than one-on-one hunting, 
animals may adopt those as simple strategies through 
natural selection process. To investigate the emergence 
of hunting patterns, evaluation of group hunting is im- 
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FIG. 1: (color online) Capturing and escaping pattern in 
group. Red and green convex objects are chasers and es- 
capees, respectively. The convex direction corresponds to a 
traveUng direction. 



portant viewpoint. 

First, we consider simple and feasible patterns in order 
to characterize chases and escapes. We intuitively clas- 
sify collective motions of capturing into several patterns. 
As a simple classification, we consider capturing patterns 
by focusing on a single escapee or a single chaser with 
its relation to opponents. The first row of Fig. [5] shows 
chaser's motions while many chasers pursue one escapee 
(A to C), and the second row shows escapee's motions 
while many escapees escape from one chaser (D to F). A 
situation of each pattern is explained as follows: A) A 
number of chasers follow an escapee in a one-directional 
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FIG. 2: Chasing and escaping processes characterized by pa- 
rameters q and p. The pattern A, B, and C represent motion 
of chasers around an escapee with corresponding values of q. 
In contrast, The pattern D, E, and F represent motion of 
escapees around a chaser with corresponding values of p. 



formation. B) Chasers surround escapees. Capturing of 
the escapees typically follow this pattern. C) One chaser 
drives an escapee into a group of chasers. This is a tran- 
sient behavior often observed leading to the pattern B. 
D) Escapees escape in a one-directional formation from 
a single chaser. E) Escapees scatter away from a chaser 
isotropically or are divided into small groups. This pat- 
tern frequently appears right after a chaser invades into 
an aggregate of escapees. F) One chaser runs after an 
escapee, while nearby escapees escape in different direc- 
tions of the chaser. 

Let us introduce a parameter q to distinguish the pat- 
tern A to C, which is assigned to escapees. At each time 
step, we focus on every escapee fc, and the + 1 chasers 



chasing the escapee k. Here, we index the nearest chaser 
as J = and other chasers as i = l,...,n^. For each 
escapee k, we define the parameter as 



qk 
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rik ■ rok, 



(1) 



where fok denotes a unit vector pointing the direction 
from the nearest chaser i = to the escapee fc, while rik 
are also unit vectors from the i-th chaser to the escapee 
k. Figure ISl^i) illustrates qk when the escapee k is chased 
by two chasers. By this parameter, the three patterns A, 
B, and C yield qk ^ 1,0 and —1, respectively. We also 
introduce the average of qk for Nt escapees which are 
chased by more than two chasers as 



Nt 



qk- 



(2) 



k=l 



We also distinguish the patterns D to F from the view- 
point of chasers rather than escapees by the similar ex- 
pression to q. At each time step, we focus on each chaser 
k, and escapees escaping from the chaser. A parameter 



{\)qk = rok- r^k 



{\\)Pk = rok ■ r^k 




FIG. 3: Illustrations for calculating parameters (i) qu and (ii) 
Pk in case of three players. 
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FIG. 4: Time evolution of parameter q for the numbers of 
chaser, Nc ~ 10, 30, 200. Initially 10 targets are randomly 
placed on the 100 x 100 square lattice. 



for the chaser k is defined in the same way as, 



1 

'k = ^ } ^ rik ■ 
^ Nc 



rok 



^ k=l 

where the denotes the number of escapees escaping 
from the chaser k. Figure Eljii) also illustrates pk when 
two escapees are escaping from the chaser k. By this 
parameter, the three patterns D, E, and F in Fig.[2]yield 
Pk 1,0 and —1, respectively. 

The above parameters q and p quantitatively distin- 
guish the patterns shown in Fig. [2l By these parame- 
ters, we can quantify group formations locally from snap- 
shots such as Fig. [1] From another point of view, this 
analysis corresponds to local structure with anisotropy. 
Therefore, we consider that the parameter is also useful 
to characterize the anisotropic structure and can be ap- 
plied to a variety of physical systems, such as amorphous 
materials Q. 

Next, we demonstrate an application of the present 
parameter to simple model of group hunting. Re- 
cently, a simple model of chase and escape in groups is 
proposed [ll[. Even this simple model shows intriguing 
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motions related to collective hunting behavior. Figure [T] 
actually shows snapshots of a simulation of the model. 
In the previous work, macroscopic quantities, the time 
for entire catch T and typical lifetime r of escapees, are 
observed but the dynamics is not understood in relation 
with the chasing spatial configuration in a group. We 
investigate details of collective motions for chase and es- 
cape in a group by q and p. 

Let us briefly explain the model in Ref. [T2|. Ini- 
tially, two types of players, named chasers and targets 
(escapees) , are placed randomly over the two-dimensional 
square lattice. A periodic boundary condition is imposed 
on the lattice. A chaser moves by one lattice unit to- 
ward its nearest target. On the other hand, a target 
tries to move away from its nearest chaser by one lattice 
unit. Here, the nearest player means the one located in 
the shortest Euclidean distance, and if there are multiple 
nearest chasers (targets), the target (chaser) chooses one 
of them with equal probabilities. Each player chooses the 
next hopping site in the following ways. When a chaser 
and its nearest target are on the same axis, the chaser 
chooses the nearest site toward the opponent, but the 
target randomly chooses one of three neighboring sites 
to increase the distance. In the other situation, chasers 
and targets choose one of two possible nearest sites with 
an equal probability in order to move closer to, or away 
from their opponents, respectively. We include exclusion 
volume effect such that players remain in the same site 
if the next hopping site is occupied. When a chaser and 
a target are placed next to each other, then the chaser 
moves to the position of the target to remove it. This 
catching rule leads to monotonic decrease of the number 
of targets. After the catch, the chaser pursues one of the 
remaining targets in the same manner. The simulation 
starts with Nj, targets and Nc chasers, and end when all 
the targets are caught. 

Figure U] shows time evolution of q for different num- 
bers of chasers with Nj, = 10. When the number of 
chasers is much larger than initial number of targets 
{Nc = 200), almost all targets are initially surrounded 
by chasers. Thus the initial q is close to 0. Since most 
of targets are immediately caught in few time steps by 
the pattern B, q remains almost cnstant. On the other 
hand, when the number of chaser is as small as the ini- 
tial number of target {Nc ~ 10), q initially fluctuates, 
and eventually approaches 1. This result indicates that 
the remaining targets are generally chased by a group 
of chasers as in the pattern A. This can be explained 
as follows. The catching event results in aggregation of 
chasers as in the pattern B. After catch, the gathering 
chasers tend to pursue the same nearest-target, leading 
to the pattern A, as the remaining targets become small. 
For an intermediate number of chaser {Nc = 30), behav- 
ior of q depends on an initial configuration. In particular, 
the initial position of players influences appearance of the 
pattern A. A fraction of samples in which the pattern A 
appears increases as the number of chaser decreases. 

Now let us turn our attention to the relation between 



time evolution of q and the number of target, NT{t). 
Figure [5] shows the time evolution of q and NT{t) for 
N^ = 10 and Nc = 5. We find that q fluctuates around 
the value 1 in most of time, but spike-like dips are ob- 
served right before catch events which correspond to the 
decrement timings of targets. This indicates that chasers 
surround targets to capture them in the patterns B and 
C. At the capturing motion, the chasers aggregate, and 
after that they form a larger group chasing a single near- 
est target. This leads to the rapid rise of g to 1. 

This parameter q also explains the crossover behavior 
of T and t 0jS cL function Nc- Here, T is defined as a 
time in which all targets are caught and r is given as 
T = Tj^o*i^T{t - 1) - NT{t))/N^. As shown in Fig. 
|6l T and t show two kinds of power-law behaviors as a 
function of Nc and kinks around Nc = 50 for N!^ = 10. 
With parameter g, we can understand now that these 
crossovers come from difference of frequency in appear- 
ance of the patterns A, B and C. As explained before 
in Fig. m q stays near (the pattern B) when Nc is 
much larger than iVj., and g ~ 1 (the pattern A) appears 
frequently when Nc becomes smaller. Here, we quan- 
tify them by observing 0: the fraction of samples which 
achieve q > 0.8 at least once until capturing all targets. 
In Fig. [SI we clearly see that the crossover point coincides 
with the point cj) becomes almost zero. In other words, 
the pattern A dominates below this crossover point, while 
almost all targets are rapidly captured with the pattern 
B or C after this point. 

The parameter q quantifies how single target is pursued 
by a group of chasers. However, it is insufficient to quan- 
tify behaviors of aggregating targets. In particular, with 
the condition N^ > Nc, a few drastic decreases in Nt 
are observed as shown later. This suggests that an aggre- 
gate of targets is caught at once surrounded by chasers. 
In such situation, the parameter p, together with q, helps 
us to understand the drastic decrease of the number of 
the targets. As an example, Fig. [7] shows time evolution 
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FIG. 5: Time evolution of q and NT{t)/N^ with A^^ = 10 
and Nc = 5. 
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FIG. 6: Dependence of T, r and (p on Nc- 4> denotes the ratio 
of samples which achieve q > 0.8. 
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FIG. 7: Time evolution of q,p and NT{t)/Nt in the vicinity 
of a drastic decay of the numb of the targets with = 50 
and Nc = 10. 

of and Nrit) for = 50 and Nc = 10. The dras- 
tic decrease of Nxit) occurs around t = 2700. Before 
the catching event, p keeps around 2 for a certain pe- 
riod. This indicates that aggregate of targets is caught 
by surrounding chasers. The upper figures in Fig.[7]show 
snapshots before and after the event. After the dras- 
tic decay, p rapidly decrease to 1 or less exhibiting the 
pattern E or F. Simultaneous occurrences of such events 
of p and drastic decay of Nxit) explains the collective 
catching motion by the surrounding targets. 

So far, we have shown results of rather small size sys- 
tem where p and q approximately represent the individ- 
ual chasing and escaping processes. In order to inves- 
tigate the average (dominant) behavior of players in a 
given condition, we investigate trend of group motion 
in group chase and escape with a large amount of play- 
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FIG. 8: Time evolution of q,p and NT{t)/N^ with N^ = 
2^^,Nc = 2",2^^2^^ The system size is 2048 x 2048. 



ers. For this purpose, we expand the system size keep- 
ing the number density of players and effectively take 
the ensemble averages of the parameters. Here, the sys- 
tem size is set to 2048 x 2048 square lattice. We note 
that individual motions of local aggregates are no longer 
captured by the parameters, however they reflect the en- 
tire spectrum. As shown in Fig. [HI (j) becomes zero for 
Nc/ L"^ > 0.01, where L denotes a linear system size. On 
the other hand, for higher density of targets, (p shows a 
finite value. Figure [8] shows time evolutions of q and p 
with N^ = 2^^,Nc = 2", 2^5 and 2^^. Even in the case 
of Nj. > Nc, the pattern B is initially dominant for the 
chased targets as g ~ 0. In addition, p ^ indicates 
that dominant escaping is in the pattern E, which means 
that most of targets can survive in a while. In the middle 
stage, plateau region appears in q and p for all Nc- In 
that region, q increase Nc decreases. In contrast, p is 
greater than but far less than 1, and targets escape al- 
most in the pattern E and occasionally with the pattern 
D. This behavior of q and p tell us that, in the plateau re- 
gion, a fair number of targets escape while chasers chase 
"decoy" targets. 

After the plateau region, q and p start to fluctuate 
as the remaining targets become small. It is notewor- 
thy that the value of q becomes closer to 1 even though 
Nc/L'^ > 0.01. This is not observed when the initial den- 
sity of targets is low. It indicates that chasers and targets 
are likely to take the pattern A because chasers chase "de- 
coy" targets at the initial and the plateau stages. The 
appearance of the pattern A contributes to relatively long 
lifetime of target as in Fig. [51 This fact corresponds to our 
previous result [l2j, which showed the lifetime of targets 
is longer in a high-density case than that in a low-density 
case for N^. 

In conclusion, we investigate group behavior of chas- 
ing and escaping by two parameters we proposed, and 
classify the motions into three patterns from a target's 
and a chaser's viewpoints. By increasing the number of 
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chasers, the kinks of power-law behaviors for the entire 
catch T and typical lifetime r of targets appear because 
of the one-directional formation, the pattern A. On the 
other hand, collective behavior of targets is also charac- 
terized. By observing time evolutions of q and p, we can 
detect the drastic catch-at-once of targets' aggregate by 
the surrounding chasers. In addition, wc find that the 
one-directional formation causes longer lifetime of tar- 
gets. 

The basic idea of the present parameters is to quan- 
tify spatial distribution of related objects from a view 
point of fine subject, and to characterize local structures 



with anisotropy. Our parameters are defined without any 
dynamical properties of objects. Then we expect that 
they can be applied to various systems. In particular, 
they have advantages to investigate structures of multi- 
component systems such as the chasing and escaping enti- 
ties, in which a characteristic property appears relatively 
locally but not in a global scale. A randomly packed me- 
dia can be another example and analyzing the local struc- 
tures could be promising by applying the parameters 14 . 
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